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Abstract 

Juggler's exclusion process describes a system of particles on the 
positive integers where particles drift down to zero at unit speed. After 
a particle hits zero, it jumps into a randomly chosen unoccupied site. 
We model the system as a set-valued Markov process and show that 
the process is ergodic if the family of jump height distributions is uni- 
formly integrable. In a special case where the particles jump according 
to a set-avoiding memoryless distribution, the process reaches its equi- 
librium in finite nonrandom time, and the equilibrium distribution can 
be represented as a Gibbs measure conforming to a linear gravitational 
potential. 

Keywords: exclusion process; juggling pattern; set-valued Markov 
process; ergodicity; positive recurrence; set-avoiding memoryless distri- 
bution; noncolliding random walk; Gibbs measure; maximum entropy 



1 Introduction 

Juggler's exclusion process (JEP) describes a system of n indistinguishable 
particles on the positive integers Z+ = {0,1,2,... }, where particles drift 
down to zero at unit speed, and after a particle hits zero, it jumps into a 
randomly chosen unoccupied site. Denote by Sn the countable space of n- 
element subsets of "Lj^. A JEP is defined as a discrete-time Markov process 
in Sn with transition probability matrix 



P{A,B) = { 



1, if O^yl, 5 = A-1, 

if G A, = - 1) U {y}, (1.1) 
0, else. 



where j4 — 1 is the set obtained by shifting the elements of A down by 
one position. A* = A \ {min ^4} is the set obtained by deleting the smallest 
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element of ^, and v = {vb) is a family of probability measures on TLj^ indexed 
by i3 G Sn-\^ satisfying 

^ v^^y) = for all B G Sn~\. (1.2) 

In general, a set A G S„ describes the heights of the particles in the system. 
The probability measure vb describes the height of a particle after a jump 
which occurred while the other n — 1 particles drifted down from config- 
uration B + 1 into configuration B. Assumption (1.2) guarantees that all 
particles jump into unoccupied sites, and therefore the number of particles 
in the system remains fixed at each transition. 

The name JEP is inspired by viewing the particles as balls that a jug- 
gler throws into a discrete set of heights according to a random pattern. In 
this context condition (1.2) amounts to a so-called siteswap juggling pattern 
where only one ball is allowed to be at hand (zero height) at any time in- 
stant. Such patterns have been studied mostly in the deterministic periodic 
setting, see e.g. Knutson, Lam and Speyer [■>] and references therein. In 
addition, Warrington [10] has computed the equilibrium distribution of a 
random siteswap pattern where throw heights are bounded and uniformly 
distributed. 

JEPs may be naturally encountered in various application areas such as 
quantum physics and engineering. For example, if particle heights are viewed 
as energy levels, a JEP can be viewed as a system of n randomly excited 
particles subject to the condition that no two particles may share the same 
energy level. Alternatively, a JEP may represent the residual job completion 
times in a manufacturing system of n machines running in parallel, which is 
operated under the condition that only one job may complete at any time 
instant. More generally, the particle heights in a JEP may be viewed as 
residual lifetimes in a discrete-time point process which may be regarded as 
a superposition of n coupled renewal processes. 

A probability measure tt on Sn is an equilibrium distribution of a JEP 
if and only if it satisfies the balance equation ttP = vr, or equivalently, 

7r{B) = tt{B + 1) + vr((A + l)U{0})z.A(5\^) (1.3) 

agb 

for all B G Sn- The following example illustrates the solution of (1.3) for 
one-particle systems. 

Example 1.1 (One-particle JEP). For n = 1 the family (i^b) reduces to a 
single probability measure 1^0 . In this case the height of the unique particle 
in the system's configuration is a discrete-time Markov process on Z_|_ with 
transition probabilities 

1, ii X > 0, y = X — 1, 
P{x,y) = { V(D{y), if a; = 0, y > 0, 
0, else. 
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This process may be identified as the residual hfetime of a discrete-time 
renewal process with interevent distribution f0 (e.g. Asmussen [i, Section 
1.2]). A direct computation shows that this JEP has an equilibrium distri- 
bution if and only if m = ^x:>qXI'q{x) < oo, in which case the equilibrium 
is unique and given by tt{x) = m~^f0([x, oo)). We note that tt = j/0 in the 
special case where f0 is geometric, a well-known fact in renewal theory. 

In this paper we will show that a wide class of JEPs have a unique equi- 
librium distribution (Theorem 2.2). However, finding an analytical formula 
for the equilibrium by solving (1.3) for a general JEP with n > 2 particles 
appears difficult if not impossible. We saw in Example 1.1 that for a geomet- 
ric jump height distribution the equilibrium is also geometric. This suggests 
that a closed-form analytical formula for many-particle systems might be 
available for some special jump height distributions. The problem is that 
condition (1.2) rules out geometric height distributions in the many-particle 
case. To get around this problem, we introduce the notion of a memo- 
ryless distribution on a subset of integers, and show that when particles 
jump according to such distribution, the equilibrium can be expressed as a 
Gibbs measure conforming to a linear gravitational potential (Theorems 3.6 
and 3.8). Such JEPs are shown to converge to equilibrium ultrafast: in a 
finite nonrandom time (Theorem 3.7). 

The model studied in this paper may be seen as a special instance of 
an exclusion process, although most literature on exclusion processes is fo- 
cused on continuous-time models with infinitely many particles, especially 
the asymmetric simple exclusion process (ASEP) on the doubly infinite inte- 
ger lattice (see Grimmett [■', Sec. 10.4] or Liggett [G, Sec. VIII] for a general 
overview, and Martin and Schmidt [ ] for discrete-time ASEP models). To 
avoid confusion, we emphasize that as an exclusion process, a JEP is very 
special because the movement of particles is deterministic, and the only 
source of randomness are the jumps from the boundary. 

The rest of this paper is organized as follows. Section 2 presents a general 
theorem for the existence of a unique equilibrium under natural integrability 
conditions on the jump height distributions. Section 3 provides a detailed 
analysis of JEPs generated by memory less jump height distributions. Sec- 
tion 4 briefly discusses two models closely related to memoryless JEPs, and 
Section 5 concludes the paper. 

2 Existence of an equilibrium 

The goal of this section is to study when the probability distribution of a 
JEP converges to a unique equilibrium distribution, regardless of its initial 
state. To rule out Markov processes with multiple recurrent classes and 
periodic dynamics, the following result gives a simple condition. 
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Lemma 2.1. Any JEP generated by a family {vb) such that 



UBiminiB")) > for all B € S, 



n 



(2.1) 



is aperiodic and has a unique recurrent class. 

Proof. Denote by G = [0, n — 1] the ground state of 5„ where the particles 
are located as low as possible. Let us write A ^ B if P*{A,B) > for 
some t > 1, and define G = {A G 5„ : G — > ^4} as the set of configurations 
accessible from G. Using (2.1) it is not hard to verify that A ^ G for any 
A € Sn, which implies that G is the unique recurrent class of the process. 
Assumption (2.1) further yields P{G, G) > 0, which implies that the process 
is aperiodic. □ 

Not all aperiodic JEPs converge to an equilibrium. For example, if each 
particle jumps to a height twice the maximum height of the other particles, 
then such a process evidently drifts to infinity. To rule out such nonphysical 
examples, we need to impose some conditions on the jump height distri- 
butions. A family of probability measures {i'b) on Z+ is called uniformly 
integrable if 



The following result shows that a wide class of juggler's exclusion pro- 
cesses is stochastically stable in the sense that as time tends to infinity, the 
probability distribution of the particle configuration converges to a unique 
equilibrium. 

Theorem 2.2. Any JEP generated by a uniformly integrable family of dis- 
tributions (ub) satisfying (2.1) has a unique equilibrium distribution vr, and 
the process started from any initial configuration converges in total variation 
to TT as time tends to infinity. 

Proof. Define the function V : Sn ^ IR+ by V{A) = max^, and denote 



where P{A,B) are the transition probabilities in (1.1). If ^ A, then 




x>K 



as K ^ oo. 



(2.2) 



PV{A)= V{A)P{A,B) 



PV{A) - V{A) = V{A - 1) - V{A) = -1. 



(2.3) 



Assume next that (z A. Then 



PViA) = ^ V{{A* - 1) U {y})uA*.i{y) 
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and because V{{A* - 1) U {y}) = mayi{V{A) - l,y), it follows that 
PV{A)= ^ (v{A)-l),,A*.i{y)+ yvA*-i{y) 

y<V{A) y>V(A) 

< v{A) - 1 + ^sup Yl vMy)- 

y>V{A) 



B£S„-i 



The uniform integrability of (i^b) implies that the last term on the right side 
above is less than ^ whenever V{A) is large enough. After combining this 
observation with (2.3), we conclude that there exists a number K > such 
that 

PV{A) -V{A) < for all A G such that y( A) > A". (2.4) 

The uniform integrability of (i^b) also implies that PV{A) < oo for all 
A £ Sn- Because the set {A £ Sn '■ V{A) < K} is finite, we conclude that 

sup PV{A) < oo. (2.5) 
AeS„:V(A)<K 

Further, Lemma 2.1 implies that the JEP is aperiodic and has a unique 
recurrent class G. As a consequence, the process is ^-recurrent (Meyn and 
Tweedie [8, Sec. 4.2]), where the measure tp on the power set of Sn is defined 
by ip{A) = \G n A\. In light of the Foster-Lyapunov drift bounds (2.4) 
and (2.5), the claim now follows by applying [n, Theorem 14.0.1] (with / 
identically one). □ 

3 JEP with memoryless height distributions 

3.1 Memoryless distribution on a subset of positive integers 

This section is devoted to analyzing JEPs generated by height distributions 

My) = { n ^ R ^^-^^ 

[ 0, y e B, 

where hsiy) = | [0, y — 1] \ -B| is the number of points in strictly less than 
y, and a £ (0, 1). The probability measure defined by (3.1) shall be called 
the memoryless distribution on with parameter a. This nomenclature 
is motivated by the following property, which states that if denotes a 
random variable with distribution ub, then the conditional distribution of 
given that > s is the same as the distribution of the random variable 

S + (,B~s- 

Proposition 3.1. For all finite B C Z^-, the probability measure vb defined 
by (3.1) satisfies 

UBitjOo) = UB[s,oo)i'B^s[t — s,(X)) for all s < t. (3-2) 
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Proof. By applying the identity hsis + t) = hsis) + hB~s{t), one can 
check by a direct computation that I'sit, oo) = a^^^^^ and I'B-sit — s, oo) = 
(yhB{t)-hB{s) £qj. positive integers s and t such that s <t. □ 

Observe that for = 0, formula (3.2) is the standard memoryless prop- 
erty on Z+, and the memoryless distribution i/0 on = is the geometric 
distribution on with success probability a. 

The memoryless distribution on can be represented in terms of a 
geometrically distributed random variable by using the following function. 
We define for a given finite A C the A-avoiding shift by the formula 

(9a(x) = min{n > : 1 [0,n] \ > x+ 1}, x G Z+. (3.3) 

The A-avoiding shift is the natural bijection from Z_(. onto the complement 
of A, and Oa{x) equals the {x + l)-th element of A'^ (see Figure 1). Some 
basic technical facts about this map that are needed in the sequel are listed 
in Appendix A.l. 

1 2 3 4 5 

\ — ^ 



eA{o) eA{i) 0a(2) 

Figure 1: The action of the A-avoiding shift for A = {0, 2, 3}. 

Lemma 3.2. Let B he a finite subset ofZ^, and let ^ be a geometric random 
variable on Z_|_ with success probability a € (0, 1). Then the random variable 
^-b(0 has the memoryless distribution (3.1) on B^. 

Proof. Because 9b maps Z+ onto B'^, it is clear that V[Ob{C) € i?) = 0. The 
claim thus follows immediately after noting that 6~^^{x) = hB{x) for x ^ B^ 
by Lemma A. 3. □ 

3.2 Noncolliding union of geometric random variables 

The noncolliding union of a finite (possibly empty) set A C Z4. and (possibly 
nondistinct) points xi, . . . , x„, G Z_|_ is denoted by 

U{A,xi,...,Xn) (3.4) 

and defined using the A-avoiding shift (3.3) recursively by 

U{A,xi) = Au{eA{xi)}, 
[/(A,xi, . . . ,Xfc+i) = U{A,xi, . . . ,Xfc) U {9u(A,xi,...,Xk){xk+i)] ■ 
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The noncolliding union of points xi, . . . , x„ E is defined by substituting 
A = into (3.4), and denoted by 

U{xi,. . . ,Xn). (3.5) 

By construction, U{A, set with [j4| + n elements which con- 

tains A but not necessarily xi, . . . , x^- Elementary properties of the noncol- 
liding union are listed in Appendix A. 2. 

Fix a parameter a G (0,1), and let ^i,^2;--- be independent random 
variables all having the geometric distribution on Z+ with success probabil- 
ity a. We will denote the n-element random set obtained as the noncolliding 
union of ^i, . . . by 

G„ = t/(6,...,en). (3.6) 

In addition, the noncolliding union of a finite B C Z+ and ^i, . . . , ^„ will be 
denoted by 

Gn{B) = U{B,ii,...,in). (3.7) 

The following two lemmas extend well-known facts for geometric random 
variables to the random sets Gn and Gn{B). 

Lemma 3.3. For any finite (possibly empty) B C such that ^ B, the 
random set Gn{B) defined by (3.7) satisfies 

E{/(G„(i?)) I i Gn{B)} = E f{Gn{B - 1) + 1) 

for all bounded functions f on S\B\+n- 

Proof. By Lemma A. 4, ^ Gn{B) if and only if min{^i, . . . , > 0. There- 
fore, 

E{/(G„(i?)) I i Gn{B)} = E {f{U{B, 6, • • • , ?«)) | 6 > 0, . . . > 0} . 

Because the geometric random variable ^fc, conditioned on being strictly 
positive, has the same distribution as + 1, and because ^i, . . . are 
independent, it follows that 

E{/(G„(i3)) I ^ Gn{B)} = Ef{UiB,Ci + i, . . . + i)), 
so the claim follows by Lemma A. 5. □ 

Lemma 3.4. For any finite (possibly empty) B C 1^+ such that ^ B, the 
random set Gn{B) defined by (3.7) satisfies 

E {f{Gn{BY) I G Gn{B)] = E f{Gn-l{B - 1) + 1) 

for all bounded functions f on SiBi+n-i- 
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Proof. With the interpretation Gq{B — 1) = B — 1, the claim makes sense 
also for n = 1. In this case the claim follows by observing that Gi{B)* = B 
on the event that = 0. To proceed by induction, assume that the claim 
holds for some n. By definition, Gn+i{B) = ?7(G„(i3), so we can split 

the event {0 € Gn+i} into disjoint events 

ni = {0(^GniB), Cn+i = 0} and = {0 G 

Let us first analyze the event Qi. On this event, Gn+i{B)* = Gn{B), so 
by Theorem 3.3 and the independence of Gn{B) and ^n+i> we see that 

E{/(G„+i(i?)*) 1 = E{/(G„(i?)) I i Gnm .3 
= E/(G„(i3-l) + l). 

Let us next analyze the event In light of Lemma A. 6, we see that 
on this event, 

Gn+l{By = U{Gn{B),Cn+ir = C/(G„(S)*, ^n+l + 1). 

Because Gn{B)* and in+i are independent, we obtain by conditioning on 
^n+i, applying the induction assumption to the map A i-^ /(C/(y4, + 
and then using Lemma A. 5 that 

E{/(G„+i(i?)*) I ^2} = E{/([/(G„(i?)*,en+i + 1)) I G Gn{B)} 
= E/(C/(G„_i(i3-l) + l,e„+i + l)) 
= E/(C/(G„_i(i?-l),en+i) + l). 

Note that 

U{Gn-l{B - l),en+l) =st U{Gn-i{B - l),^™) = - 1), 

because ^„ and ^n+i have the same distribution and both random variables 
are independent of Gn-i{B — l). (We use the notation X =st Y to mean that 
the random elements X and Y have the same distribution.) This implies 
that 

E{/(G„+i(i?)*) I J^2} = E f{Gn{B - 1) + 1). (3.9) 

By combining (3.8) and (3.9), we conclude that the claim is true for n + 1. 
This completes the induction step, and thus also the proof. □ 



3.3 One-step evolution 

Assuming that the height distributions are memoryless as in (3.1), we see 
with the help of Lemma 3.2 that the JEP may be constructed using the 
noncolliding union defined by (3.4) according to 

f - 1, ^ Xt, 

Xt+i = { , , ' (3.10) 

I c/(x;-i,6+i), OG^i, ^ ' 
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where ^i,^2)--- are independent geometric random variables in TLj^ with 
success probabihty a, independent of the initial state Xq. 

The following key technical result describes the one-step evolution of the 
JEP with memoryless height distributions, when started at a random initial 
state distributed according to GkiB) defined in (3.7). 

Lemma 3.5. If Xq =st Gk{B), then the JEP with memoryless height dis- 
tributions at time 1 is distributed according to 

^ ^ f Gk{B-l), O^B, 
' ''\Gk+,{B*-l), GB. 

Proof. Assume first that ^ B. Fix a bounded function /. Then by first 
using (3.f0) and then Lemma 3.3, we find that 

E {/(Xi) I ^ Xo} = E {fiGkiB) - 1) I ^ GkiB)} 
= Ef{Gk{B-l)). 

Further, because in (3.f0) the random variables ^t+i and — 1 are inde- 
pendent, 

E {/(Xi) I G Xo} = E {f{U{Gk{By - 1, a+i)) I G Gk{B)} , 

so by Lemma 3.4 and independence, 

E{/(Xi) |0 G Xo} = Ef{U{Gk-i{B - l),Cfc+i)) 
= Ef{Gk{B-l)). 

By combining these two observations, we conclude that the claim holds when 
O^B. 

Let us next analyze the case where G Then with the help of 
Lemma A. 6 we see that 

Xi =st U{Gk{By 

= [/(c/(i?*,ei + i,...,^fc + i)-i,efc+i) 

= [/(c/(s*-i,6,...,4),Cfc+i) 

= Gk+i{B*-l). 

□ 

3.4 Ultrafast convergence to equilibrium 

By applying the key result. Lemma 3.5, we may identify the equilibrium of 
the memoryless JEP. 

Theorem 3.6. The JEP generated by the memoryless height distributions (3.1) 
has a unique equilibrium which has the same distribution as the noncollid- 
ing union Gn = f7(Ci, . . . ,^n) of n independent geometric random variables 
defined in (3.6). 
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Proof. Assume that the JEP is started at a random initial state distributed 
according to G„. Then by Lemma 3.5 apphed with = 0, it fohows that 
Xi =st Gn- By the Markov property, it follows that Xt =st Gn for all t. The 
uniqueness of the equilibrium distribution follows from Theorem 2.2. □ 

Lemma 3.5 also allows to make a stronger conclusion on the precise evo- 
lution of the JEP started at an arbitrary nonrandom initial state. Consider 
a JEP started at a nonrandom initial state A with n elements. We will 
denote by the set obtained by deleting k smallest elements from A, and 
by Tfc the time instant when the A;-th smallest particle from the initial con- 
figuration has just jumped, so that ti = min^ + 1 and = min^fc_i + 1 for 
k = 2, . . . ,n. The following result shows that the memory less JEP reaches 
its equilibrium at the nonrandom time instant when all particles have 
jumped once. 

Theorem 3.7. Let {Xt)te'i,^ be a JEP generated by memoryless height dis- 
tributions (3.1) and started at a nonrandom initial state A ^ Sn- The 
distribution Xt at any time t can be represented by 



Xt — St < 



A — t, for t < Ti, 
U{Ak -t,S,i,. . . ,^k), if Tk<t < Tfc+i for some 1 < /c < n - 1, 
U{ii,. . . ,in), fort>Tn. 



Proof. First, Xt = A — t for t < ri by the definition of ti. When t = ti, we 
see using the representation (3.10) that Xt = U{A* — t,^t) =st U{Ai —t,S,i). 
Thereafter, by repeatedly applying Lemma 3.5, we have Xt = U{Ai — t,^i) 
for all Ti < t < T2. Note that € — t for t = T2 — 1. Therefore, again 
by Lemma 3.5, Xt =st U{A2 — t, ^i,'^2) for t = T2. By continuing the same 
way, we conclude that 

Xt =st U{Ak - t, 6) • • • )^fc) 

for all t such that < t < t^+i for some 1 < k < n — 1. Finally, at the time 
instant t = Tn — I, Xt =st U{An-i —t,(,i,... ,^„_i), where the set An-i — t 
only contains the zero element. By one more iteration of Lemma 3.5, we 
conclude that Xr„ =st U{^i, . . . ,Cn)- This distribution remains invariant by 
Theorem 3.6. □ 



3.5 Distribution of the equilibrium 

Having found a stochastic representation for the equilibrium of the JEP with 
memoryless height distributions as the noncolliding union of n independent 
geometric random variables, it is natural to ask whether there exists a nice 
formula for the distribution of this random set. The following result provides 
an affirmative answer. 
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Theorem 3.8. The random set Gn = U{£,i, . . . defined by (3.6) is dis- 
tributed according to the Gibbs measure 

^lp{B) = Z-le-P"^''\ BeSn, (3.11) 

where f3 = — logo, the potential energy of B is given by H[B) = X^^g^ 
and the normalizing constant (a.k.a. partition function) is given by 

n a 

ei^ 



= n rr- (3.12) 



C' 
k=l 

Proof. The claim fohows by substituting A = ^ into the more general The- 
orem 3.9 below. □ 

Consequently, the memoryless JEP attains maximum entropy among all 
n-element random sets with bounded expected sample mean (e.g. Cover 
and Thomas [2, Sec. 12.1]). As a byproduct, Theorem 3.8 also provides a 
computationally efficient way to generate samples from distribution (3.11). 
A direct computation using (3.11) shows that for the memoryless JEP in 
equilibrium, the expected sample mean height of the particles is given by 

the probability of the ground state by 

n 

F{Gn = [0,n- 1]) = "Q 

k=l 

and the probability that zero is occupied by 

P(0 eGn) = l- a". 

By Birkhoff 's ergodic theorem, the last formula is also equal to the long-run 
average jump rate in the system — a simple way to prove the formula is to 
observe that the minimal element of Gn is geometrically distributed with 
parameter a" due to Lemma A. 4. 

We will next state and prove a generalization of Theorem 3.8, which 
allows to compute the nonequilbrium distributions of a memoryless JEP 
started at an arbitrary nonrandom initial state, as described in Theorem 3.7. 
Denote the family of supersets of a finite (possibly empty) set A with + n 
elements by 

SA,n = {B CZ+: Bd A, \B\A\=n}. 

Theorem 3.9. The random set Gn{A) = U{A,^i, . . . ,^„) defined by (3.7) 
is distributed according to the Gibbs measure 

fiAAB) = Z-^ e-'^^^(^), B G SA,n, (3.13) 
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where 13 = — logo, is the constant defined by (3.12), and Ha{B) = 

"^x^BXA^A^x) is the ^-neglecting potential energy of B, where 

hAix) = \[0,x-l]\A\ (3.14) 

is the number of elements in strictly less than x. 

Proof. Consider first the case n = l. \{ B = AiJ {x} for some x € A^, 
then Gi{A) = B \i and only if ^a('^i) = x. The probability of this event is 
(1 — a)a^'^^^^ by Lemma 3.2. Because Ha{B) = h^ix) for such B, the claim 
follows for n = 1. 

Assume next that the claim is true for some n > 1. Fix a set B D A 
such that \B \ A\ = n + I. Then 

P{Gn+i{A) = B) = PiGn{A) = B\{x}, 9B\{,}iCn+l) = x) 

xeB\A 

= ^ F{Gn{A)=B\{x})P{6B\{x}{Cn+l)=x). 
x&B\A 

Because Ha{B \ {x}) = Ha{B) — hA{x) for all x G A^, it follows by the 
induction assumption that 

V{Gn{A) = B \ {x}) = a^^(^)-''^(-). 

Further, 

V{eB\{x}{in+i) = x) = (l-a)a'^™(-) 
by Lemma 3.2. Note also that hB\{x}{x) = hsix) and 

hA{x) = hB{x) + \[0,x-l]n{B\A)\. 

By combining these observations, we find that 

P(G„+i(^) = i?) = (l-a)Z-Ja^^(^) (3.15) 

xGB\A 

To compute the last sum, observe that | [0, x — 1] fl (i? \ vl) | equals k — 1 when 
X is the k-th smallest element in B\A. Therefore, 

^ ^-|to,x-i]n(B\A)| ^ ^'a-(fc-i) = ^ (3.16) 

x€B\A k=l " 

By combining (3.15) and (3.16), we conclude that the claim is true for 
n + 1. □ 
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4 Related models 



4.1 Bounded JEP with uniform jumps 

Consider the n-particle JEP where jump heights are bounded by a constant 
M > n, and the particles jump into unoccupied sites in [0, M — 1] uniformly 
randomly. In this case Warrington has shown that the equilibrium 
distribution is given by 

^,{B) = Z-},l[{l + \[x + l,M-l]\B\), (4.1) 

xeB 

where the normalization constant Zn^M = { m^i^^ } is a Stirling number of 
the second kind (Stanley [: , Corollary 2.4.2]). This equilibrium may also be 
expressed as the Gibbs measure 

where /? = 1 and 

Hm{B) = - J] log (1 + |[x + 1, M - 1] \ 

Recall that the memoryless JEP reaches its equilibrium when the initially 
highest particle jumps for the first time. The same is not true for the 
bounded JEP with uniform jumps. To see this in the one-particle case, 
observe that for n = 1 the equilibrium distribution (4.1) reduces to 

/x({x}) = Z-Ij{M -x), X G [0, M - 1], 

whereas just after each jump, the particle is uniformly distributed in [0, M — 
!]• 

4.2 Finite ASEP with a reflecting boundary 

Consider a continuous-time asymmetric simple exclusion process on Z_|_ con- 
sisting of n identical particles, where each particle attempts a jump x i— )• x+1 
at rate A and a jump x i— > x — 1 at rate rj, independently of the other par- 
ticles (e.g. Grimmett [3, Sec. 10.4]). An attempted jump takes place if the 
targeted site is unoccupied and belongs to Z^; otherwise the attempt is 
suppressed. 

Assume that < A < 77, and let Y be the Markov jump process in 
the Weyl chamber Wn = {x £ Z" : xi < X2 < ■ ■ ■ < x„} which keeps 
track of the particle locations. If Y denotes the corresponding free process 
in Z" where particles are allowed to jump on top of each other, then the 
components of Y are independent M/M/1 queues, and so the equilibrium 
distribution of y is a product of n independent geometric distributions on 
Z+ with success probability X/rj. Further, because Y is reversible, we obtain 
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the equilibrium of Y be truncating the equiUbrium of Y to the Weyl chamber 
Wn (e.g. Kehy [4, Cor. 1.10]). As a consequence, we see that the equiUbrium 
of the corresponding set- valued process in Sn is the Gibbs measure (3.11) 
with (3 = — log(A/?7). 

Thus, the memory less JEP and the finite ASEP with a reflecting bound- 
ary have the same equilibrium distribution. 

5 Conclusions 

This paper introduced a juggler's exclusion process (JEP) as a special discrete- 
time exclusion process on the positive integers, where the motion of particles 
is deterministic outside the boundary, and particles hitting the boundary 
jump into random unoccupied sites. We showed that a JEP is ergodic if the 
jump heights are uniformly integrable and satisfy a natural aperiodicity con- 
dition. The main part of the analysis was devoted to JEPs where the jump 
distributions satisfy a set-avoiding memoryless property. Such JEPs were 
shown to converge to equilibrium in finite nonrandom time, where the equi- 
librium is a Gibbs measure corresponding to a linear gravitational potential. 
The proof of this result was based on representing the time-dependent sys- 
tem configuration in terms of independent geometric random variables. As a 
byproduct, the proof yielded a fast computational way to generate samples 
from such Gibbs measures. 

There are two open problems related to this model that we find par- 
ticularly interesting. First, do we obtain an ergodic JEP if we assume that 
jump distributions are bounded in but not uniformly integrable? Second, 
is it possible to write down an analytical formula for the equilibrium of a 
general ergodic JEP using techniques from renewal theory, as in the spe- 
cial case of one-particle JEP? The model analyzed in this paper also leaves 
room for many generalizations. For example, one might consider analogous 
models with a continuous time parameter and fully random particle motion 
as in ASEP models. Another interesting direction is to look for continuous- 
time processes where particles move in continuum and jump into random 
locations not too close from the other particles. 
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A Technical details 
A.l Set-avoiding shift 

Recall that Oa is the ^-avoiding shift defined by (3.3) in Section 3.1. 
Lemma A.l. For all finite A C Z+ and all x € 

eA+i{x + i) = eA{x) + i. 

Proof. Because the equations 

min{n > : | [0, n] \ yl | > x + 1} 

= min{n > : | [0, n + 1] \ A | > x + 1} + 1 

and 

|[0,n + l]\(^ + l)| = |[0,n]\Al + l 
hold for all A C and all x € Z+, we obtain 

eA+i{x + l) =min{n>0: | [0, n] \ (A + 1) | > x + 2} 

= min{n > : | [0, n + 1] \ (^4 + 1) [ > x + 2} + 1 
= min{n > : | [0, n] \ A | > x + 1} + 1 
= OAix) + 1. 

□ 

Lemma A. 2. For any finite nonempty A C 'Lj^ and any x E 

{X, X < minA, 
9a*{x + 1), x>minA. 

Proof. The first claim is that 

min{n > : | [0, n] \ A | > x + 1} = x 

whenever [0, x] \ A = [0, x]. Because 

I [0, n] \ ^ I > X + 1 = I [0, x] I = I [0, x] \ A I 

implies that n > x, the claim follows. 
The second claim is that 

min{n > : 1 [0, n] \ ^ | > x + 1} = min{n > : [ [0, n] \ ^* | > x + 2} 

whenever x > minA. Note that | [0,y] \ A | = | [0, y] \ j4* | — 1 whenever 
y > minj4. If 

y e {n>0:\[0,n]\A\>x + l}, 

then I [0, y] \ A I > X + 1 = I [0, x] I > [ [0, x] \ ^ [, so y > X > minA. By the 
note above, \ [0,y]\ A \ = \ [0,y]\ A* \ - 1, i.e. 

y e {n > : I [0,n] \ A* I > x + 2}. 

The other direction is similar, and the sets to be minimized are the same. □ 



15 



Lemma A. 3. For any finite A C "L^ and any x G A^, 

9j^{x) = hA{x), 

where hA{x) = \ [0, x — 1] \ A\. 

Proof. Now the claim is that 

min{n > : [[0,n] \ A\ > \[0,x - 1]\A\ + 1} = x 

whenever x ^ A. In this case | [0, x] \ A | = j [0, x — 1] \ j4 | + 1, so the claim 
reduces to 

min{n > : | [0, n] \ A | > \ [0,x]\A\} = x, 
which is trivial. □ 

A. 2 Noncolliding union 

In this section we summarize basic technical properties of the noncolliding 
union map U{A,xi, . . . defined by (3.4) in Section 3.2. 

Lemma A. 4. For all finite A C Z+ and all xi, . . . ,Xn G 

min U{A, xi, . . . , x„) = min(yl U {xi, . . . , Xn})- 

Proof. Let us first prove the claim for n = 1. We may assume that A is 
nonempty, because otherwise the claim is trivial. If xi < min(^), then 
dAixi) = xi, which implies that m.mU{A, xi) = xi. If xi > min(^), then 
dA^xi) > min(^), which implies that mmU{A,xi) = min(A). Therefore, 
we conclude that the claim holds for n = 1. The general case follows by 
induction, due to the recursive construction of U{A, xi, . . . , Xn). □ 

Lemma A. 5. For all finite A C Z+ and all xi, . . . ,Xn G 

U{A + l,Xi + l,...,Xn + l) = U{A, Xi,...,Xn) + l. 

Proof. Note first that by Lemma A.l, 

u{B + 1, y + 1) = + 1) u {eB+iiy + 1)} 

= {B + l)u{eB{y) + l} 
= U{B,y) + l 

for all finite B C and all y € Z+. Therefore, the claim is true whenever 
n = 1. If the claim holds for some n > 1, then by the induction assumption 

U{A + l,Xi + l,...,Xn + l) = U{U {A + l,Xi+l,...,Xn + 1), X„+i + 1) 

= U{B + l,Xn+l+l), 

where B = U{A, xi, . . . , Xn). The claim now follows by applying the prop- 
erty with n = 1. □ 
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Lemma A. 6. For any finite nonempty A C and any x S 

A, X < mill A, 
U{A*,x + l), x>min^. 

Proof. First, if x < miiiA, then = x by Lemma A. 2. Therefore, 

U{A,x)* = (A U {x})* = A. Next, if x > minA, then 6'a(x) = Oa'^x + 1) 
by Lemma A. 2, and minC/(A,x) = minA. It fohows that 

C/(A, x)* = A* U {9a{x)} = a* U {Oa* (x + 1)} = UiA*,x + 1). 

□ 
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